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Abstract 

In this paper we show that the performance of the quantum adiabatic algorithm is determined 
by phase transitions in underlying problem in the presence of transverse magnetic field T. We 
show that the quantum version of random Satisfiability problem with 3 bits in a clause (3-SAT) 
has a first-order quantum phase transition. We analyze the phase diagram 7 = 7(r) where 7 is an 
average number of clauses per binary variable in 3-SAT. The results are obtained in a closed form 
assuming replica symmetry and neglecting time correlations at small values of the transverse field 
r. In the limit of F = the value of 7(0) ~ 5.18 corresponds to that given by the replica symmetric 
treatment of a classical random 3-SAT problem. We demonstrate the qualitative similarity between 
classical and quantum versions of this problem. 
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I. INTRODUCTION 



From the early years of computing it was realized that some problems are inherently 
intractable. This intuition has been quantified by the theory of computational complexity 
proposed by Cook Most problems of practical interest can be roughly divided into two 
classes: P and NP-complete. Problems in the former class can be solved on a computer in 
time that scales only polynomially with the size of an instance of the problem. Solution to 
NP-complete problems can be verified in polynomial time, but it is believed that it cannot be 
found in polynomial time on a classical computer. Solving NP-complete problem typically 
requires exponential time which makes then intractable. All NP-complete problems are 
equivalent; if it were possible to solve one NP-complete problem in polynomial time, it 
would be possible to apply the same algorithm to solve all NP-complete problems. 

It is not known yet whether NP-complete problems can be solved efficiently on a quantum 
computer. Shor's algorithm for a quantum computer can solve in polynomial time the 
number factoring problem that is presumably hard for a classical computer. It is not, 
however, in the NP-complete class (and its decision variant is in P). 

Just as research in classical computing focus on worst-case complexity has been super- 
seded by the analysis of typical case complexity and application of general-purpose algo- 
rithms like simulated annealing, similar changes take place in the field of quantum com- 
puting. A substantial interest has been generated by a general-purpose quantum adiabatic 
algorithm (QAA) proposed by Farhi and coworkers [sl, y. 

In its simplest form QAA is applied to problems where underlying variables can have 
only two values, and a solution is given by a N-bit binary string. It corresponds to a time- 
dependent Hamiltonian that slowly changes from a simple form (for which the ground state 
can be constructed exactly) to a complex form that describes an instance of NP-complete 
problem. The initial Hamiltonian is typically chosen to correspond to a uniform magnetic 
field r applied along x direction 

TV 

i^in = -r5^af, (1) 

i=l 

SO that the ground state in {(xf } basis is a symmetric superposition of 2^ binary vectors 

^i^ = ^ Yl 1*1^2 ■■■stv). (2) 

{s,=±l} 



The final Hamiltonian is chosen to be diagonal in {o"f } basis with diagonal elements cor- 
responding to the cost function values for particular assignments of binary variables {si} 



At intermediate times the Hamiltonian is a linear combination of initial and final Hamilto- 
nians. We choose to write this in the following form 



In the beginning of the algorithm F — oo so that the second term dominates and a ground 
state has a simple form. At the end of the algorithm F = and the ground state corresponds 
to a solution of the NP-complete problem that minimizes cost function. If F is lowered 
sufficiently slowly, adiabatic theorem tells us that the system will remain in its ground state 
with a high probability. The second term is referred to as a driver term because its presence 
in otherwise diagonal Hamiltonian is to enable spin flips. The algorithm is quite similar to 
the simulated annealing (SA) algorithm in this respect. In QAA transverse field F replaces 
temperature T. The dynamics of these two problems are quite different. 

The fundamental limitation of the SA algorithm is the critical slowing down at the point 
of classical phase transition. The system may become trapped in one of the local minima 
surrounded by high barriers and that requires exponentially long time to escape. For the 
QAA the metric of the algorithm performance can be given in terms of the eigenstates ^^^(F) 
and eigenvalues -Efc(F) of the interpolating Hamiltonian H(r). The rate at which F(t) can be 
lowered obeys the strong inequality ;i|f (t)|(^o| (iif/(iF|^o) < (-^i(r) --Eo(F))2. A sufficient 
condition for the success of QAA requires a polynomial scaling with the problem size of a 
minimum gap Ai^min = niinre(o,oo) {Ei{T) — Eq(T)) between the two lowest energy levels of 
H{r). The danger is that there are points of avoided crossing ((see Fig.Q)) at which the gap 
can be exponentially small. At such point the branches of Eq^T) correspond to wavefunctions 
with exponentially small overlap. This effect has been demonstrated analytically in a toy 
problem P, 0] as well as for the case of NP-complete problem Positive- 1-in-K-sat with large 
number of bits in a clause, ^ 1 . 

An equivalent measure of the algorithmic complexity is the second derivative of the 
ground state energy d^£'o/dF^. This quantity can become exponentially large at the point 




(3) 



H = ^ E[{si}]\si ■ ■■sn){si ■■■sn\ - Py^gf. 



(4) 
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FIG. 1: Two lowest eigenvalues of the Hamiltonian H ^ vs the interpolatmg parameter r = 
1/(1 + r). In simulations i^gn corresponds to a randomly generated instance of the version of 
Satisfiability problem Exact Cover 4] with the number of bits A^=10. Insert: different plots show 
the T-dependencies of the gap between the two lowest eigenvalues of H for various random instances 
of Exact Cover with A^=10. Dashed lines indicate the approximate asymptotical position of the 
avoided-crossing point in the limit N ^ oo. 

of avoided crossing where it is oc 1/Ai?min. If the minimum gap AE'min shrinks to zero in the 
thermodynamic hmit N ^ oo for some F = then the ground state energy Eq^T) possesses 
a cusp (cf. Fig. H)). This cusp signals a first-order quantum phase transition. Indeed, the 
ground state energy is precisely the free energy in the limit /5 = 1/T — oo: 



EoiT)= lim F(/3,F) 



hm -lnTr|e-^^(r)|. 



(5) 



Time evolution has to be slowest at the point of the first-order quantum phase transition, 
reminiscent of a critical slowing down at the point of classical phase transition. 

We will be working with an ensemble of random instances of NP-complete problem and 
recast it as essentially a spin glass model. Although quantum phase transitions in infinitely- 



connected models of spin glasses have been extensively studied j8|, |9|, |10|, [ll|, there are 
no results to our knowledge on dilute quantum spin glasses. At the same time most of 
NP-complete problems are closely related to dilute spin glass models. 

The central quantity we will be computing is the disorder-averaged (instance-averaged) 
free energy {F{f3)). We will classify the quantum phase transition as random first-order 
or random second-order as deduced from the non-analyticity in {F) vs T. In classical case. 



random first-order phase transition is typically associated with exponentially hard problems. 
It is therefore tempting to suggest that the onset of the first-order quantum phase transition 
is indicative of the exponential, or stretched-exponential scaling law of the runtime 

of QAA with determined by the distribution of the tunnelling times between the valleys 
of an energy landscape. 

Within current formalism we are limited to this qualitative picture. The determination 
of the minimum gap requires more elaborate methods. One would naively expect that finite- 
size scaling analysis can be used to derive the scaling of the minimum gap with A^. Observe 
though that, if disorder is relevant and assuming the first-order phase transition, the position 
of the minimum gap fluctuates with disorder around some value Tc (see Fig. Q). If we fix 
a value of T and perform the disorder-average, we will have {Ei{T) — Eq{T)) > 0{l/\^). 
Minimized over F, it is still 0(1/ \/N) even though true minimum gap can be arbitrarily 
small. 

The paper is organized as follows. In Sec. |n] we introduce the random K-SAT problem 
- the model that we study in this paper. In Sec. IIIII we describe exact results from the 
random graph theory and empirical results on satisfiability transition. K-SAT problem has 
been extensively studied in the classical limit for zero temperature. This section is devoted 
to these results as well as effects of finite temperature. The Sec. IIVI formulates the K-SAT 
problem in quantum case in the presence of transverse magnetic field. Subsequently two 
approximations are made to make the problem tractable: the so-called replica-symmetric 
ansatz that assumes absence of long-range correlations is made, and the static ansatz that 
assumes that any correlations time correlations are due to the transverse magnetic field 
alone. The expression for the free energy and the self-consistency equation for the order 
parameter are derived. In the Sec. we consider the limit of small transverse magnetic field 
F and derive the simphfied self-consistency equations for this case. Sec. EH is Conclusion. 



II. K-SAT PROBLEM 



As a test case for QAA algorithm we consider the satisfiability problem that happens to 
be the first problem to be associated with NP-complete class l|. More precisely, we will 
work with a variant of satisfiability problem - the K-SAT that places a constraint on the 
number of variables that can appear in a clause. The K-SAT is known to be NP-complete 



5 



for ^ 3. This means that any NP-complete problem can be reformulated as an instance 
of 3-SAT. We will concentrate on this particular case {K = 3). The benefits of working 
with K-SAT is that the problem Hamiltonian is local (K-local to be exact, but can easily be 



recast as a 2-local Hamiltonian; see 1^ for application to QAA algorithm); and also that 
random K-SAT is intrinsically related to random hypergraphs and is amenable to methods 
of statistical physics. 

An instance of K-SAT is a boolean formula in conjunctive normal form (CNF), i.e. a set 
of clauses 

= Ci A Cs A • ■ ■ A Cm (6) 

where each clause represents a logical OR of K literals, each literal representing either one 
of the variables xi,. . . ,xn or a its logical NOT. Below is an example of possible K-SAT 
clause for = 3: 

a;2 V ^5 V xu (7) 

(where x represents logical NOT of x). 

The formula is said to be satisfiable if and only if there exists an assignment of boolean 
variables {xi} such that all M clauses comprising a formula are satisfied (evaluate to true) 
at the same time. Presented with an instance of K-SAT formula, algorithm must determine 
if it is satisfiable and if it is, find the appropriate assignment of variables. 

We follow the standard recipe for trivial mapping of constraint satisfaction problems 
to problems in statistical mechanics by defining the energy function to be proportional to 
the number of violated clauses. The constant of proportionality is chosen to be 2. The 
motivation for this choice will become evident later. The energy can be written as follows 

E = 2Y,^iJmlS^^^nJm2S^^,) " " " ^ ( Jmi^S.^,- ) • (8) 

TO 

where we replaced boolean variables with spins Sj = ±1, the value s = — 1 corresponding 
to true. 6{x) is the Heaviside function 6{x) = 1 if x > and 6{x) = if x ^ 0. Variables 
imi, . . . ., imK are indices of variables that appear in a clause, and J^i, . . . , JmK = ±1 describe 
whether the literal that appears in a clause is negated (J = —1). 

We will work with an ensemble of random instances of K-SAT with fixed number of 
variables A^, number of clauses M and number variables per clause K. This merely means 
that indices imp are independent and drawn uniformly at random from the set 1, 2, . . . , and 
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variables Jmp are independent random variables that are equal to +1 or —1 with probability 
50 % . 

The formula is satisfiable if the ground state energy equals zero, and it is unsatisfiable if 
it is greater than zero. 

The following phenomenon has been discovered Q|. If we fix the ratio of clauses to 
variables 7 = M/N, the energy is almost surely for 7 < 7c in the limit N 00 and 
E > for 7 > 7c. The value of 7c is independent of in the limit of large A^, apart from 
small corrections. Note that the ground state remains exponentially degenerate across the 
transition. The exact value of 7c is an interesting problem in combinatorics Q|; rigorous 
determination of 7c has not yet been accomplished. 

The threshold phenomena in random graphs were known since pioneering work of Erdos 
and Renyi 17[. Random K-SAT problem has intimate connections to the theory of random 
graphs. One can picture an instance of random K-SAT as a random hypergraph. Variables 
correspond to vertices of the hypergraph and clauses correspond to hyperedges that join K 
vertices. One threshold phenomenon that was examined in 12| in the context of random 
graphs but can be extended to random hypergraphs, is the appearance of giant component. 
In physics it is known as percolation transition. When 7 < 1/K{K—1) a random hypergraph 
is with high probability is a collection of disjoint clusters of typical size 0(1), with the largest 
cluster not exceeding O(logA^) in size. All clusters (with few (i.e. 0(1)) exceptions) have 
a form of trees, that is contain no loops. For 7 > 1/K{K — 1) the same picture holds, but 
there appears a giant (of size 0(A^)) supercluster. The structure of the giant supercluster is 
very complex. It contains a large number of loops, but the length of these loops is relatively 
large: 0(log A^). It is in the giant supercluster that the complexity of the problem is buried. 
Indeed, small isolated clusters can be efficiently tackled using divide-and-conquer approach. 
Moreover, since they are all tree-like, the assignment of variables that satisfies corresponding 
clauses always exists and can be found in linear time. 

Although phase transition is evident in geometric properties, owing to long length of the 
loops in the supercluster, no abnormal behavior in the space of solutions is seen immediately 
following the phase transition. The formula remains satisfiable well past the percolation 
transition, nor any correlations between spins in the supercluster that are far away from each 
other are detected. One must note that such correlations do appear at so called dynamic 
transition 7^, at which point finding the solution becomes difficult (for simulated annealing 
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FIG. 2: Example of hypergraph corresponding to random instance of 3-SAT. Vertices represent 
variables and traingles represent clauses involving K = 3 variables. 

algorithm) . 



A i log T 
N 




FIG. 3: Complexity of the instance as the function of 7 = M/N . Problem becomes exponentially 
hard for 7 > 7^ and the complexity peaks at satisfiability threshold 7c. 

The performance of algorithm is subject to similar threshold phenomenon. Typically, 
algorithms require just a polynomial (often linear) time for 7 less than some critical value 7^ 
and require an exponential time for 7 > 7^. One can introduce the normalized complexity 
limjv^oo log^) where T is the time it takes to solve a problem (units of T are unimportant 
in ^ 00 limit since the logarithm is taken). It is a self-averaging quantity and hence 
a function of 7 alone (in the usual sense, with high probability a random formula has this 
complexity). This function equals zero for 7 < 7^ and for 7 > 7,^ it is non-zero. The exact 
value of 7(i and the shape of the complexity function is not universal, but rather algorithm- 
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dependent. Simulations suggests that for most common algorithm the complexity peaks 
at 7c that is the point SAT/UNSAT transition. Indeed, the satisfiability threshold is the 
place where one intuitively expects to find the hardest (and most interesting) problems. The 
exact determination of the complexity function is invaluable as a non-empirical means for 
comparing the asymptotic efficiency of various algorithm. This is particularly useful in the 
field of quantum computing, where empirical research is impossible since no prototype of 
quantum computer exist and its simulation on classical computer is not feasible for large 
instances of the problem. An important benchmark of the algorithm is 7^ itself (larger is 
better) as it marks a region where the problem can be solved very efficiently. An alluring 
(but not necessarily impossible) goal is designing an algorithm that has 7d ^ 7c- 

Although the value of •jd is algorithm-dependent, a particular value that we denote 7rsb 
is universal in a sense that it should be the same for all local search algorithms. It corre- 
sponds to the point where the energy landscapes qualitatively changes. There appears an 
exponentially large number of local minima, i.e. a set of states corresponding to E = be- 
comes disconnected. In contrast for 7 < 7rsb a set of states with E = is connected. This 
transition can be seen in the behavior of the free energy. This is the point where so-called 
replica-symmetric ansatz breaks down. This will be elaborated on in the next section. 



III. REVIEW OF CLASSICAL RESULT 



In a number of articles 18|, |l9|, |20|, |2l|, |22| satisfiability transition has been studied in K- 
SAT problem using replica method. Statistical properties of physical system are completely 
determined by its free energy 



-T\nZ 



-Tln^e" 

{Si} 



(9) 



If the Hamiltonian depends on the disorder, we must also average the free energy over 
disorder configurations. Note that it is not appropriate to attempt to compute the disorder 
average of Z. Indeed Z is exponentially large and small o{N) fiuctuations of E amplify in 
Z = e~^^ . Therefore 



-/3(F) 



(10) 



Note that although the free energy E is 0{N), fiuctuations due to disorder are only 0{yN) 
and can be neglected in the thermodynamic limit. Disorder averaging of E is just a useful 



trick, since to the leading order in N, same value should be obtained for F for almost all 
possible disorder realizations (all but a fraction that goes to with increasing A^). 

Replica method accomplishes the averaging as follows, n non-interacting copies of the 
system are prepared, all having the identical disorder configuration. This is indicated by 
attaching an additional replica index a = 1, ... ,n to each variable s,; new variables are 
labeled s°. The partition function of such system equals Z", where Z is the partition function 
of non-replicated system. Disorder average of the replicated system is then performed, which 
is easy to accomplish since summation over all possible spin configurations {sf} can be 
done after the disorder averaging, which eliminates the difficulty of computing the partition 
function that explicitly depends on the disorder. Once analytical expression for (2'") is 
obtained, the disorder-averaged free energy (F) = — T(lnZ) is computed via analytical 
continuation using the identity 



The Hamiltonian for the K-SAT problem has been introduce in the previous section (Eq. 
|H)). In the limit of zero temperature satisfiable phase is characterized by (E) = 0; in the 
unsatisfiable (UNSAT) phase (E) > 0. 

In the simplest approximation it is assumed that the symmetry of the Hamiltonian with 
respect to the interchange of replica indices a is not spontaneously broken. The physical 
interpretation of this is the absence of long-range correlations. For randomly chosen variables 
Si and Sj, replica symmetric ansatz implies 



(randomly chosen sites are with high probability at least O(logA^) away from each other; 
correlations are absent in the limit A^ oo). 

Owing to identity (fT^ statistical properties (correlations) in the thermodynamic limit are 
completely determined by specifying average magnetizations of individual spins = (sj). 
The order parameter of the system is the histogram of magnetizations of spins P{m). In 
practice it is convenient to define effective fields hi = Ttanh~^mj and choose the histogram 
of effective fields as an order parameter. An expression for the free energy can be written 



{lnZ)=4-{Zn 



(11) 



n=0 



(SiSj) = {Si){sj) 



(12) 
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out in terms of P{h). 



F = \ f d/iln(2cosh/5/i) / — P(cj)(lnP(c^) - 1) 

1-e-^ 



(13) 



^ j=l,2,3 

Varying this expression with respects to P{h) yields a self-consistency equation of P{h) 
which can be solved numerically. This equation is essentially an iterative procedure for 
determining correct magnetic fields for clusters without loops (trees) combined with the 
Poisson distribution for the number of branches in a random tree. 

For high connectivities 7 and low temperatures T solution to the self-consistency equation 
yields two solutions. A solution that maximizes the free energy should be chosen. Note that 
this is in contrast to standard of choosing a solution with the smaller free energy for first- 
order phase transition in pure systems. This reversal is standard feature of replica method 
and the rationale is discussed in [2^. 

For zero temperature the solution is drastically simplified. P{h) has the form of a series 
of delta-function peaks at integer values of h 

P{h) = Y,Pk5{h-k). (14) 

k 

These magnetic fields can be interpreted as follows. Non-zero values oih{\h\ ^ 1) correspond 
to frozen variables. Indeed hi ^ 1 corresponds to nii = +1 and hi ^ —1 corresponds to 
rrii = —1. The appearance of finite fraction of such frozen spins, or backbone, signals the 
beginning of unsatisfiable phase. The absolute value of hi {\hi\) then indicates the increase 
in the number of violated clauses if the frozen spin is flipped. 

Since the problem is symmetric with respect to sign flip Si —>■ —Si, P{h) is necessarily 
symmetric. Introduce variable q 

(l = ^Pk (15) 

k=l 

Then obviously po = I — 2q since all p^s must add up to 1. The self-consistency equation 
can be written in the following form 

1 - 2g = e-=^^'''/o(37g2) (16) 
Pk = e-'^^"h{3^q') (17) 
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Equivalently, these equations can be derived as follows j2J]. Let the fraction of frozen 
variables be 2g, with half of these being polarized to +1 and another half to —1. Next, we 
add + 1-st (cavity) spin to the system together with extra clauses. The number of extra 
clauses is Poisson-distributed with parameter 87. Each extra clause involves a cavity spin 
and two randomly chosen variables of A^-spin system. The clause forces a certain value for 
the cavity spin only if both variables other than a cavity spin are frozen to a value that 
equals corresponding Jmi- The probability of this occurrence is g^. Every such clause gives 
a contribution u = J mi to the effective field h of cavity spin. This contribution is +1 or 
— 1 with probability of 50 % each. If the probability that the clause contributes u = +1 (or 
u = —1) is the number of clauses attached to cavity spin that contribute u = +1 (or 
u = — 1) is Poisson with parameter |7g^. Therefore, we can write the probability that the 
effective field of the cavity spin is k: 



The right-hand side is evaluated to be e"^'*''^ /fc(37g^). When the effective field equals zero, 
the cavity spin is not frozen. Since the properties of A^-spin system and + 1-spin system 
should not differ in the thermodynamic limit, the probability of this is 1 — 2q. Thus we 
obtain the self-consistency equation on q. 



This equation does not determine q uniquely; when both trivial solution (g = 0) and a 
non-trivial {q 7^ 0) are present, the correct value of q is chosen by examining the expression 
(fT^ for the free energy and choosing a solution that maximizes the free energy. Although a 
non-trivial solution appears at 7 ^ 4.67, it does not become stable until 7^ ~ 5.18 

At finite but small temperature this picture is modified as follows. A series of delta- 
function peaks in P{h) are broadened and acquire a finite 0{T) width (similar broadening 
occurs in quantum case for small F, see figure IH)). The values of q and integrated probability 
weights around integer values of h remain the same. Effective fields for spins that are not 
frozen h ^ acquire 0{T) corrections so that the magnetization m = tanh j3h = 0(1) and 
spins in the backbone {\h\ > 1) acquire 0{T) correction because of clauses connected to the 
backbone (see figure [71). 




'l—m,k 



(18) 



l,m 



l-2g = e-3^'''jo(37g'). 



(19) 
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Examination of the free energy also shows that the value of 'jc{T) increases with increasing 
temperature. The following phase diagram is obtained. 




FIG. 4: Simple phase diagram obtained within RS ansatz 

Unfortunately, the predicted value of T = UNSAT transition of 7c ~ .... is far from 
the experimental value of 7c ~ 4.2. The reason for the discrepancy is that the replica- 
symmetric ansatz breaks down well before the satisfiability transition. Using an improved 
1-step replica-symmetry breaking ansatz (1-step RSB) the break-down of replica symm etry 
is put at 7(i ~ 3.92 and the satisfiability transition is computed to occur at 7c ~ 4.27 jl9. 22 1. 

The central assumption of RSB ansatz is that for 7 > 7d the ground state (set of all sat- 
isfying assignments) is broken in an exponentially large number of isolated islands separated 
by large barriers. The identity (siSj) = {si){sj) no longer holds unless the thermal average 
is restricted to a particular local minimum (or pure state in statistical physics lingo). Rather 
than being described by individual magnetizations m^, the system ought to be described by a 
histogram of magnetizations taken over all possible pure states a. At the level if 1-step RSB, 
the following identity holds in the thermodynamic limit. For sites i and j chosen at random 
(and as a result, typically infinitely far away from each other) P{mimj) = P{mi)P{mj), 
where P{mi) is a histogram of magnetizations at site i, and P[mi,mj) is a histogram of 
pairs of magnetizations at sites i and j. 

Note that this ansatz is broken when once one goes beyond single step of replica symmetry 
breaking. However, it is widely believed to be exact up to 7 = 7c; thereafter, additional 
steps of RSB are required. 

Due to the complex nature of the ground state in RSB phase, simulated annealing al- 
gorithm should take an exponentially long time to converge to a solution (although other 
specially designed algorithms can still work efficiently for 7 > 7^;). A correct phase diagram 
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should take the following form. The exponential slowing down that affects the performance 
of simulated annealing algorithm occurs at RS/RSB boundary. 

At 



PM 




''RSp,----' 
/ UNSAT 



d c 



Y 



FIG. 5: Correct phase diagram obtained by considering effects of replica symmetry breaking (RSB). 
RSB phase separate "easy" and unsatisfiable phases. 

IV. QUANTUM 3-SAT 

For quantum adiabatic evolution algorithm, the original Hamiltonian is rewritten with 
operator af replacing classical spins and a driver term — FS^af is added 



We expect level crossings and a critical slowing down to occur at phase transition indicated 
by non-analyticities of quantum partition function 



Note that we will let /3 — > oo and examine the behavior of Z as a function of transverse field 
F. 

To evaluate the partition function we use the Trotter lattice. Equation ()21|) can be 
rewritten as 




(20) 



Z = Tre-'^^ 



(21) 



(22) 



We can use Baker's identity 



= eVe-5[^'^l+- 



(23) 
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where additional terms involve higher-order commutators. In our case both A and 13 are 
0(1/M) hence the last term in (|^^ is e*^*^^/^^^^ and does not contribute to (|^^ in the limit 
M ^ oo. 



Tr 



-{f3/M)H,i[{a-}] (I3r/M) 



M 



(24) 



To evaluate this expression for each of M factors we insert basis states corresponding to 
all 2^ projections of a-. We denote basis states by vector Si^t where t = 0, . . . , M labels 
a factor in which the basis states appear. Since we are taking a trace we are imposing a 
periodicity condition = Si^m- Since Hc\ is obviously diagonal in this basis, first term is 
simply written as 

and using e"°^'' = cha + a^sha, second term can be written as 

l^ch— j ginth(/3r/M)E,{i-*».t^..t+i)/2_ (26) 

The partition function of system with quantum spins is thus written effectively as a 
partition function of classical system with NM classical spins. 

Z = (chTAt)^*^ ^ e-'^tT.tHcil{sr,t}]+inth{rAt)j:i^{i-s,,tSi,t+i)/2^ (^27) 

where we have denoted At = j3/M for readability. We are obviously interested in the limit 
M ^ 00, /3 ^ 00, At ^ 0. 

In what following we shall focus on the case with 3-bits in a clause (3-SAT). As we have 
described earlier, the classical Hamiltonian Hd for 3-SAT takes the following form: 

H4{s,,t}] =2J2^iJalS^MJa2SJMJa3SkJ (28) 

a 

Variables ia,ja and ka are drawn randomly out of 1, . . . , A^; variables Jai, Ja2, Jas are equal 
to or —1 with probability of 1/2. A particular realization of these variables is referred 
to as disorder. 

In replica method we prepare n identical copies of the system so that the partition function 
becomes Z". In practice, spin variables are augmented with a replica index a = 1, . . . ,n 

= (chTAt)"^^ e"^*^".'^^«^(-^''i''-)^('^'^2^^a)e(.^a3SftJ+inth(rAt)E<..,,t(i-<t<t+i)/2 (^29) 
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Note that after disorder average is taken, Hamiltonian becomes symmetric with respect to 
permutations of sites labeled by index i (but not to permutation of replica indices a or 
Trotter indices t). Due to this permutation symmetry, mean field theory should be exact. 
We introduce the following variables 

^(K}) = ^En^K*'<]- (30) 

i a,t 

This counts the fraction of sites for which spin assignment equals certain vector cr". By 
definition S{o-j"}c({cr"}) = 1. Disorder-averaged summand of (j^^ can be expressed entirely 
in terms of c({crf }). The summand can be written as exp(Vi + V2), where the first term is 
due to clauses 

V,[c{ct)]=^N Y1 c(p)c(o-)c(T)(e-2^*S-'^(^i''*)^(-^^'^?)^(-^«-*")> (31) 

the remaining average being only over the signs of Ji, J2, J3. 
Second term is due to the transverse magnetic field 

VM^)] = iVlnth(rAt) J2 E ^^4^ (^^^ 

{af} a,t 

In the limit N 00 sum over all possible realizations of c(cr) can be replaced by a continuous 
integral. The entropic term (due to multiple realizations of {s^t} with identical c(cr)) takes 
a simple form 

5[c(<t)] = -iV^c(o-) In c(<t). (33) 

The partition function can be written as 

(Z") = (chTAt)"^^-^ J [dc(o-)]e^[^('")]+^^[^('")]+^=^[^('")l (34) 

Since terms in the exponential are 0{N), the value of the integral is determined by its 
saddle-point. 



A. Replica-Symmetric Ansatz 



A major simplification is obtained if we assume that c({(j"}) is symmetric under permu- 
tation of replica indices a. For simplicity, the following substitution is made 

'{rt} 



CUT 



[dH{T)]P[H{r)]\{ ^ (35) 



16 



where the integral is over all possible single-site Hamiltonians (i.e. all possible vectors of 
size 2^ corresponding to all spin configurations on Trotter lattice). 

This has the desired symmetry property and is more amenable to taking n ^ limit, 
since c((t) is now encoded by the probability distribution P[H{t)]. Maximization over all 
possible c(cr) is replaced by maximization over all possible P[H{t)]. 

The entropic term is computed as follows 



-^c(^)lnc(^) = -A^[e(^)p 



(36) 

p=l 



Substituting replica-symmetric ansatz for c{cr) we obtain 



- y ^ - - ^ - - ^ ^^\Z[H,{ct)]...Z[H,{ct)] 

(37) 

where we have denoted Z[H{a)] = ^{^t} 

Taking the limit n ^ we use the fact ~ 1 -|- n In x and keep only contribution linear 
in n: 



J][c(c7)f = n /"[di/i(o-)]P[i/i(<7)]...[di/j,(<7)]P[i/j,(o-)]lnZ[i/i(<7) + ... + i/^(o-)038) 

-np J[dH{(T)]P[H{(T)]lnZ[H{(T)] (39) 

Note that H{a) is defined only up to a constant. Without losing generality we can require 
that X^{o-t} ^{^) ~ 0- With this constraint, we define the Fourier transform of P[H[cr)]: 

P[uj{cr)] = j [dii'(<7)]e-'^{<^*}'^('^)^('^) (40) 

Owing to the constraint X^i^t} ^i*^) ~ 0' Fourier transform P[a;(<T)] is invariant under shift 
io{cr) uj{(t) uoq. With the aid of this Fourier transform the first term (which is basically 
the integral of convolution of P[H[(t)]) can be rewritten as 

J[dH{a)]\nZ[H{(T)] J [duj{(T)]P[uj{(T)]Pe'^i'^t}'^^''W^^ (41) 

(Note: a factor of (27r)~^*^ is implicit in [da;(<T)]) 

This makes differentiation over p trivial. The final result is 

S^-N j [dH{cT)] In Z[H{a)] j [du;{(T)]P[u;{fT)]{\nP[u;{(T)] - l)e'^{^t}^(''W'r} ^43) 
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The remaining contributions are obtained straightforwardly 



1=1,2,3 



n p[HMmM], 

1=1,2,3 

(43) 



where the averaging is done over the signs of Ji, J2, J3 and is defined as follows 

Z;[H^{a),H2{a),H;{a)] = ^ eHi(p)+i/.(,.)+//3(x)-2AtE.e(Jip0e(J2<x0e(J3n)_ (44) 

{pt},Wt},{rt} 

And for V2 we obtain the following expression 



V2 = N J [dH{a)]P[H{a)] 



Mlnch(rAt) + ln^e^(")+^°*^(^^*)^*(^-"'"*+^)/2 - lnZ[i7(cr)] 

Wt} 

(45)- 



B. Static Approximation 

Since working with the most general form of H{a) is intractable, we make an ansatz 

H{(t) = hAt ^ (Tt + In th(rAt) (rtcxt+i)/2 (46) 

t t 

with a single parameter h. This describes an isolated spin subject to time-independent 
external magnetic field h. Alas, any dynamic effects are neglected within this approximation. 
This is similar in spirit to the static approximation made in solving infinitely-connected 
model in jsj. Neglecting dynamic effects still permitted to obtain a qualitative picture of 
the spin glass phase. It is widely believed that static approximation works best in the limit 
of small r. We specifically consider that limit in the next section. 

The important simplification of present approach is that the order parameter becomes 
a simple function P{h). Moreover, we can now take a limit M ^ 00 further simplifying 
calculations. To the lowest order in At 

d/iln[2ch/3Vr2 + /i2] / — e^^^P(cu)[lnP(cj) - 1] (47) 

-00 J—Qo 

-N\n{TAt) f dhP{h) (3Tth ((3VT^ + hA , (48) 

J —00 

where P{uj) is the Fourier transform of P{h). 
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For Vi we obtain the following expression 



{\uZ^{hi,h2,h^))j^ j^j^- \nZi{hi) 

1=1,2,3 



n p{hMK 

1=1,2,3 



(49) 



where and Z^ can, respectively, be expressed as follows 



Zi{h) = Tre^('^*^+^'^-\ 
Zsihi, h2, /la) = iYe^Ei=i,2,3('^<°'i^+ro-ix)-2/36»(Jio-i;j)e(J2iT2^)e(J3'73^ 

For V2 the following expression is obtained to leading order in At. 

/oo 
dhP{h)(3rthpVT^ + h^ 
-00 



(50) 
(51) 



(52) 



which exactly cancels At dependence in S. Hence, as expected, for sufficiently large M, 
solution is independent of M. 

Observe that correct classical free energy is obtained if F is set to zero for finite p. For 
purely quantum result we set P — +00, so that the free energy F — —{S + Vi + V2)/P 
becomes 



F 

N 



dhVr^ + / — e''*^'*P(a;)[lnP(u;) - 1 



{A{Jihi, J2/i2, Jshs)) 



Jl,J2,J3 



1=1,2,3 



(53) 

JI P{hi)dK (54) 



=1,2,3 



where 



A(Ji/ii, J2/i2, Jshs) = Ai, 



.1=1,2,3 1=1,2,3 1=1,2,3 



(55) 



with Amax denoting the largest eigenvalue of corresponding 8x8 matrix. The result depends 
on Ji, J2, J3 only through products J\hi, J2/12, Jshs- It is subsequently averaged over signs 

of Ji, J2, J3. 

Differentiating the free energy with respect to P{h) allows us to write self-consistency 
equation for P{h). 



P{oj) — exp37- 



/ [(A(Ji/ii, J2/12, ^3/i3))j„j„j3 - E.=i,2,3 \/r^TIFl e-''M/iiP(/i2)d/i2P(/i3)d/i3 



VF2 + K^Qi^^'dh 



(56) 
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Replacing Jihi by hi the averaging can be thrown out 
P{uj) = exp 37 



/ A'ih, h2, h) COS Ujhdh ^(^^) Y^-"^) dh2 



(57) 



+ /i^ COS ujhdh 

where A'(/ii, /i2, /is) = A(/ii, /i2, /^a) - Ei=i,2,3 + hj. 

From this form it is also evident that P{oj) is real and hence P{h) is symmetric P{h) = 
P{~h), as it should be by the symmetry of the model. 

In high-F, low-connectivity phase we expect a solution to the self-consistency equation to 
be unique, whereas in high-connectivity, low-F phase two solutions are present. To determine 
which solution to take, one must examine the free energy F and choose a solution that 
maximizes its value. By rewriting (|54|) with the aid of self-consistency equation a somewhat 
simpler expression is obtained 



F 

N 



27 f A'{h,h2,h) n Pih^)dK 
1=1,2,3 



-37 



A'(+oo, h2, /is) + A'(-oo, /t2, /is) 



n Pih,)dh, 



i=2,3 



(58) 

VF2 + h^P{h)dh 



V. SMALL r LIMIT 

We examine the solution to the self-consistency equation in the limit of small F. We 
seek a solution in the form of a series of peaks of width 0(F) centered at integer values of 
h (see Fig. EJ. This behavior is similar similar to what happens in classical case for small 
temperatures. 




+2 



+3 



FIG. 6: Probability distribution of effective fields for small transverse field F. 

We can use the degenerate perturbation theory to evaluate A{hi, h2, /is)- It will invariably 
involve finding the maximum eigenvalue of some matrix corresponding to nearly degenerate 
levels. To facilitate the discussion we denote by L2{a), L^{a,h), L4^{a,b,c), LY{a,b,c), the 
largest eigenvalues of certain matrices A^j (j = 2, 3, 4, 7) given in Appendix. 
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The reason for finite width of peaks is very similar to that for classical case. Variables in 
isolated clusters and variables in the backbone connected to free variables attain a small 0(r) 
to their effective fields (see figure [7j). Since trees of arbitrary depth are technically possible, 
exact solution requires exact diagonalization of arbitrary large matrices. The nature of the 
static approximation that we made replaces effects of distant clauses by an effective field; 
due to this the maximum size of the matrix to be diagonalized is 2^ x 2^ - the same as for 
isolated clauses. Since we are doing this procedure in self-consistent manner, it is better 
than simply truncating the expansion in size of the cluster. 

The matrices we consider above are submatrices of 8 x 8 matrix that correspond only 
to those rows and columns that involve only combinations of three spins that keep clause 
satisfied. For isolated clause this leaves only 7 combinations. In the limit of F = these 
correspond to the 7 degenerate energy levels; finite F lifts the the degeneracy and diago- 
nalization of 7 X 7 matrix is required. When one or more spins in a clause are frozen, the 
degeneracy is smaller. The above-mentioned matrices enumerate all possible degeneracies 
for K = All possible combinations of effective fields that give rise to thses expressions 
are given below. 




FIG. 7: Spins in the backbone (black area) together with some clauses connected to them. The 
effect of the clauses is OiV) corrections to effective fields of spins in the backbone. 

With the aid of the notation that we just introduced, we will be able to write the resulting 
expressions in compact form. We focus on the interaction term in the expression for the free 



21 



energy 

V^N-f I A{hi,h2,h3)- J2 + n P(h,)dh, (59) 

L i=l,2,3 J 1=1,2,3 

For the remainder of this section when we say h ^ k where k is an integer, we actually mean 
that {h-k)/r ^0(1). 

Let us evaluate the expression in parentheses A' = A{hi,h2,hs) — Z]j=i 2 3 V^'^ + ^1- 
Since the expression is symmetric we can assume hi ^ h2 ^ without losing generality. 
To first order in F the following expressions are obtained for all possible cases: 

1. hi < 0, \hi\ > F. In this case A' = 0. 

2. Neither /ii fa nor hi ^ 1. Classical expression holds in this case. A' = — 2min(l, hi). 



3. hi fti 0and/i2 ^ F. No degeneracy. Simple answer A' — —hi — ^/T^~+l^. 

4. hi ~ 0, /12 ~ Oand/13 ^ F. Triple degeneracy. For A' we obtain the following 
expression: 



A' = TL^{-hi, -h) - F 5^ 1/1 + hi (60) 

i=l,2 

where hi = hi/T. 
5. hi 0, /i2 ~ 0, /i3 0. 7-fold degeneracy. For A' we have 

A' = VL^{-hi, -h2, -h) J2 + (61) 



=1,2,3 



where hi — hi/V. 

6. hiP^l and(/i2 - 1) > F. Double degeneracy. A' = -2 + TL2{hi) - Thi. 

7. /ti a; 1, /i2 ~ 1 and(/i3 — 1) ^ F. Triple degeneracy. 

A' = -2 + VL^{hi, h) - T{hi + h), (62) 

where hi — {hi — 1)/F. 

8. /ii 1, /i2 ~ 1, ^3 ~ 1- 4-fold degeneracy. 

A' = -2 + TL^Chi, h2, /13) - r(/ii + h2 + hs), (63) 
where hi — {hi — 1)/F. 
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The self-consistency equation for P{h) can be written in the following form 



TOO / 5V/n \ iu,h^]L 
J-oo \SP{h) J ^ 



where K{uj) = J y/T^'+We^^'^dh. It is convenient to write 



K{u) = -^QiVu), (65) 



where Q{x) can be expressed using modified Bessel function of the second kind 



\x 



Q{x) = -^K,{\x\). (66) 



/27r 

Also observe that Q{0) = 1. With this replacement we can write 

'1 d^ fSV/N\ 



~ 1 /-^ 1 d^ SV/N\ 



e'^^dh. (67) 



We seek a solution in the form of a sequence of peaks around integer values oi h {h ^ k), 
each peak having a width 0{T). Write 



^c) - f E A. (^) . 

fc=-oo ^ ^ 



(68) 



Total probability weight around h ^ k can be expressed as Pk = Pk{h)dh. 
The interaction term V has the form proportional to 

j A\hi,h2,h3)P{hi)P{h2)P{h3)dhidh2dh. (69) 

Computing a variation 6V/6P{h) leaves only two probability distributions in the integral; 
say P{h2) and P{hz). 

First, we assume that neither /i nor h 1. Note that the two remaining fields in the 
integrand: /i2 and hz are necessarily integers, lih < /i2, ^3 we know that A' = 9{h) min(l, /i), 
which, differentiated twice over h gives just zero. 

II h > h2, and h — h2^T,h — h^^r, then 5V/SP{h) is actually independent of h 
and does not survive double differentiation over h either. Note that this holds even if /i pa 1 
(but either /i2 ~ or /la 0). 

Therefore, the only non-zero contributions are from h ^ and /i ~ 1. Moreover expres- 
sions of the form P{h)dh and J;^^ P{h)dh can be replaced with constants q and 
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Q ~ Pi (q — Ylk=iPk) since we have seen aheady that varying with respect to P(h) in these 
expressions leads to /i-independent terms. 

It is convenient to rewrite V in terms of Pkih)- We write 



V/N = 27g^ + 7ry' , 

where V takes more complex form 

y - I Lr{-h,-h2,-h) H Po{k)dk + 3q j L,{-h,,-h2) J] PoChi)dh, 

j=l,2,3 i=l,2 

+3q^ J hiPoChi)dhi-3{po + qf J Vl + ~h?Po{h)dh 

+ I Li{hiMM n Pi{hi)dhi + ?,{q-p{) I L^ihiM) J[ Pi{hi)dhi 

i=l,2,3 1=1,2 

+^{q-Pif j L2{h)Pi{hi)dhi - 3q^ J hPi{h)dh. 



(70) 



(71) 



We can use the identity 5V/5P{h) = 5V/5Po{h) ior h ^ and 5V/5P(h) = 5V/5Pi(h) for 
h ^ 1. We have shown that all other terms do not contribute. 

For the partial derivatives we obtain the following expression (also taking into account 
the fact that Po(h) = Po{-h). 
1 6V/N 



^^^SPoCh) 



1 6V/N 
37r 6P,{h) 



L7{kkMPo{h2)Po{h3)dh2dh3 + 2q J L^{h~h2)Po{h2)dh2 



Wh - {Po + qfL2{h) 



lS, h2, h3)Pi{h2)Pi{hs)dh2dh3 + 2{q - pi) J Ls{h, /is) A(/i2)d/i2 



+ iq-p,fL2Ch)-q^h 



We write the self-consistency equation in the following form 



P{u;) — exp 



37 
Q{u) 



FoiTu) + ^ (e-Fi(ru;) + e--F*(ru;)) 



where we have denoted 

Fo(a;) = 



-oo 
oo 



1 / 1 SV/N \ 
2dh^ \3jr 6PoCh)) 
1 / 1 6V/N\ 
2dh? V37r SPiCh)) 



e-'^^dh, 



''^^dh. 



(72) 

2 

(73) 

(74) 

(75) 
(76) 



Since for all functions L„ we have the identity \^L{h, . . .)\^_'^ — 1, we obtain Fo(0) = —q^ 
and Fi(0) = q^. 
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We now perform the Fourier transform of P{uj). Since P{uj) is modulated periodic func- 
tion with cUmax ~ l/r? its Fourier transform is necessarily a series of spikes of width F. This 
justifies our previous ansatz. Let us compute the probability density for h ^ k. 

PkCh) = r e^'^-'^^^'exp [FoiTu) + Re {e'^ F.iTuj)}] ^. (77) 

Using the smallness of F and the identity J^^ Qacosui+ikui^ _ can rewrite to leading 

order in F 

J-oo \ Q{^) J 

Note that the integrated probability weights Pk = J Pk{h)dh are obtained by substituting 
u) = in the integrand 

Pk = e-''^^'h{3^q'). (79) 

This is precisely Eq. For k = using pq = 1 — 2q a. self-consistency equation identical 
to that of F = 0, T = is obtained. Therefore, the value of q is unchanged to leading order 
in F. 

Since -Fo('^) and Fi{uj) are given solely in terms of Po{h) and Pi{h), k = and k = 1 are 
sufficient to provide a closed system of equations 

PoCh) = r e3^^o(^)/«(^)/o ( ^^^^^) e^^^ — , (80) 
J~oo V Qy^) J 27r 

P,Ch) = re^^^°(^)/Q(^)Ji P^'^^(f^' le-^^ — ■ (81) 
J-oo V Q{^) ) 27r 

For some critical value of 7 both trivial (g = 0) and a non-trivial (g 7^ 0) solutions coexist. 
While the appearance of non-trivial q and its value are not sensitive to F for small values of 
F, the point where the non-trivial solution becomes stable is. 

Also observe that since the integrand is real, all Pk{h) are symmetric Pk{h) = Pk{—h). 
While this is expected for Po{h), such symmetry for Pi{h) and others is likely only approx- 
imate, higher order contributions in F should make Pk{h) 7^ Pk{—h) for k 0. 

Determining the stability of non-trivial solution is accomplished with the aid of Eq. ()59p. 
Substituting our ansatz for P{h) we obtain 

L^^ = 37g2 - 47g3 + ^ \k\p, - 2^TV' - ^^TV" + T f Vl + h'Po{h)dk (82) 

k 
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where V has been defined in ()71|) and V" (coming from A'(+oo, . . .) term) can be written 

as 

V" = J LsChiMPiChi)PiCh2)dhidh2 

+ 2(g-pi) j L2Chi)PiChi)dh~2q j hPi{h)dh. (83) 

For r = the last three terms disappear and substituting pk = Ik{'i'l<f) "we obtain 

^ = iq' (l - e-^^''^ (/o(37g^) + M,{?>^q'))) (84) 

The value of 7 = 7co ~ 5.18 where this expression (with g 7^ determined from self- 
consistency equation) becomes positive is the point where non-trivial solution becomes sta- 
ble. We now calculate the phase transition line 7 = 7c(r) to the leading order in F (here 
7(0) = 7co). For 7 slightly larger than 7co, classical expression can be written as 



^ 0.0706. (85) 

7c (0) 



The remaining terms in F/N are linear in F. Write them in the form —TV{q): 



6) 



V{q) = 27 /" LrChMM n PoChi)dh^ + 6^q I L^{hiM) \{ Po{k)dhi 

i=l,2,3 '' i=l,2 

+27 I L^{hi,~h2,h3) II PiChi)dk + 6l(\ + q-pij j L^{hiM) W Pi{hi)dhi 

i=l,2,3 ^ ^ 1=1,2 

+67(9 -Pi) Q + g-Pi^ J Vl + h^Pi(h)dh - (1 + 67(1 - qf) j Vl + h^PoCh)dh 

Let Vo(go) denote V{q) computed at 7 = 7c(0) with Po{h), Pi{h) solving self-consistency 
equation with g 7^ 0; denote Vo(0) similarly computed V{q) at 7 = 7c for trivial solution 
g = 0. 

To leading order in F, the point 7c where nontrivial g 7^ solution becomes stable is 

7.(r) = o.(o) + rM_m. (87) 



VI. CONCLUSION 



In t his p aper we have extended the classical treatment of phase transitions in K-SAT 



n t nis p 

HQ 



I8L I20I |21[ to the quantum domain for the case of i^' = 3. Although infinitely-connected 
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quantum spin glass models have been studied, no studies of dilute spin glasses have been 
performed to the best of our knowledge. 

While infinitely-connected models have small {0{1/\/N) or smaller) couplings, dilute 
glasses are characterized by strong (0(1)) couplings which means that perturbation expan- 
sion cannot be truncated. Due to the limitations imposed by the structure of disorder, we 
have only been able to solve the problem within static approximation using replica symmetric 
ansatz. 

What we have observed is that the quantum limit F > 0, T = is qualitatively similar 
to the classical limit F = 0, T > (although quantitative results and analytical expressions 
are quite different). The order parameter in the limit of small F takes the form of series of 
peaks of width 0(F), just as in the classical case it takes the form of series of peaks of width 
0(T). At F = 0, T = we recovet the phase transition at the classical value 7c (0). For 
small but finite F the value of 7c increases linearly, 7c = 7c(0) + OF, and the expression for 
the constant O is given in the closed form in terms of the quantities computed at F = 0. 

Much of the similarities can be explained away by the fact that we ignored dynamic effects. 
Incorporatin g th e dynamic effects changes the phase diagram of Sherrington-Kirkpatrick 



model 



jratm g ti i 

Qyin 



However, static approximation is assumed to work very well in the limit 



of small F. For the problem at hand this seems to be the only feasible limit. A bigger 
concern is that we have completely ignored the effects of replica symmetry breaking (RSB). 
We expect that the location of the djTiamic transition for the QAA algorithm should be the 
same as for simulated annealing, since it is given by F = 0, T = 0. Working in the regime of 
small T or small F within RSB will enable us to compare performance of these algorithms 
for 7 ^ 7d. 

Another suggestion for future work is K-XOR-SAT problem. It can be solved on a 
classical computer in polynomial time, but becomes exponentially hard for the simulated 
annealing algorithm. Owing to simple structure of its energy landscape, many exact results 
have been obtained for this problem [2^. Quantum version of i^-XOR-SAT in the limit of 
small F might be much easier to solve. Note that many of properties of K-SAT are found in 
K-XOR-SAT, example being a single level of replica symmetry breaking in a certain range 
of 7. 

We have worked within replica symmetric formalism developed by Monasson j26[ that 
uses a functional order parameter. Alternative method of working with dilute glasses that is 
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closer in spirit to the treatment of SK model has been developed by Viana and Bray [27| and 
it incorporates a sequence of various order parameters. The bridge between this approaches 
have been developed by Kanter and Sompolinsky j^]. A more consistent approach to 
making static approximation would be to ignore time dependence in a sequence of VB-like 
order parameters along the lines of ^ and convert the result to the form that uses functional 
order parameter. We have not tried to reconcile that approach with our treatment. It is 
interesting to see if these approaches are equivalent, and if not, how the answer changes. 
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APPENDIX 



Matrices Aij introduced in the Sec. |3have the following form: 



a 1 
1 -a 



a + h I 1 
1 a-h 
1 -a + h 



( a+b+c 1 

1 a + h — c 



1 

V 1 








1 



a — h + c 




1 






\ 



-a + h + c J 



Ml 



b + c 


1 


1 





1 





1 


a + b — 


c 


1 





1 


1 





a-b^ 


-c 1 











1 


1 


a — b — c 








1 











-a + 6h 


-c 1 





1 








1 


—a + b 








1 





1 
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